A criterion for the nature of the superconducting transition 
in strongly interacting field theories : Holographic approach 
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It is beyond the present techniques based on perturbation theory to reveal the nature of phase 
transitions in strongly interacting field theories. Recently, the holographic approach has provided 
us with an effective dual description, mapping strongly coupled conformal field theories to classical 
gravity theories. Resorting to the holographic superconductor model, we propose a general crite- 
rion for the nature of the superconducting phase transition based on effective interactions between 
vortices. We find "tricritical" points in terms of the chemical potential for U(l) charges and an 
effective Ginzburg-Landau parameter, where vortices do not interact to separate the second order 
(repulsive) from the first order (attractive) transitions. We interpret the first order transition as the 
Coleman- Weinberg mechanism, arguing that it is relevant to superconducting instabilities around 
quantum criticality. 

PACS numbers: 



Interactions between vortices contain information on 
the nature of the superconducting transition. They 
change from repulsive to attractive, decreasing the 
Ginzburg-Landau parameter re, the ratio between the 
penetration depth of an electromagnetic field and the 
Cooper-pair coherence length Combined with ei- 

ther the e = 4 — d expansion or the 1/N approximation in 
the Abelian-Higgs model [4j, |5[ , one finds that the nonin- 
teracting point for vortices at re = re t (~ I/a/2) is identi- 
fied with the tricritical point, where the nature of the 
superconducting transition changes from second order 
(re > re t ) to first order (re < ret) [J]. Quantum corrections 
due to electromagnetic fluctuations are the mechanism, 
referred as the fluctuation-induced first-order transition 
@ or Coleman- Weinberg mechanism [||. 

The situation is much more complicated when cor- 
related electrons are introduced. In particular, super- 
conducting instabilities are ubiquitous in the vicinity of 
quantum critical points 7], where quantum critical nor- 
mal states are often described by strongly interacting 
conformal field theories. Although one can integrate over 
such interacting fermions, the resulting effective field the- 
ory contains a lot of singularly corrected terms for Higgs 
fields, which originate from quantum corrections due to 
abundant soft modes of particle- hole and particle-particle 
excitations near the Fermi surface [1, Q. Furthermore, 
the Fermi surface problem turns out to be out of control 
[Tol fill since not only self-energy corrections but also 
vertex corrections should be introduced self-consistently. 
It is far from reliability to evaluate effective interactions 
between vortices in this problem. 

Recently, it has been clarified that strongly coupled 
conformal field theories in d-dimension can be mapped 
into classical gravity theories on anti-de Sitter space in 
d + 1-dimension (AdS<j+i) [H, EH- This framework has 
been developed in the context of string theory, refereed as 
the AdS/CFT correspondence. See Ref. [3]' for a review. 
Immediately, it has been applied to various problems 



beyond techniques of field theories: non-perturbative 
phenomena in quantum chromodynamics (AdS/QCD or 
holographic QCD) llal . non-Fermi liquid transport near 
quantum criticality |l6l4l8l | and superconductors [l9l I20I ] 
in condensed matter physics (AdS/CMP), and etc. 

In this letter we propose a general criterion for the 
first-order superconducting transition based on the holo- 
graphic approach. We take the holographic superconduc- 
tor model [2fJ as an effective low-energy model in the dual 
description for certain classes of strongly interacting field 
theories. The asymptotic vortex solution [21] turns out 
to play a central role in the nature of the superconduct- 
ing transition. We suggest "tricritical" points in terms of 
the chemical potential for U(l) charges and an effective 
Ginzburg-Landau parameter, where vortices do not inter- 
act to separate the second order (repulsive) from the first 
order (attractive) . We interpret the first-order transition 
as the Coleman- Weinberg mechanism [1,0], arguing to be 
relevant to superconducting instability around quantum 
criticality. 

We start from the holographic superconductor model 
in AdS4 with radius L 
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where the complex scalar field is decomposed into the am- 
plitude 77 and the phase tf>, and is the bulk gauge po- 
tential with the field strength F = dA. re p is the Planck's 
constant. In this work we set e = 1 and m 2 = —2/L 2 
and consider the probe limit. The background metric is 
given by 

dz 2 \ 

* 2 f{z)dt 2 + dx\ + dx\ + —— , (2) 
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with f(z) = 
by T= |2. 
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The Hawking temperature is given 
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Equations of motion read 

V 2 r 1 -m 2 r 1 -j 1 Ql = 0, 
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rfQ v 



-L 2 V^, 



(3) 



where = — V^4> is the gauge invariant superfiuid 
four- velocity and is its field strength. is d^d^cf) 
which can be replaced with delta functions for centers of 
vortices. 

Now we calculate effective interactions between vor- 
tices. The effective interaction will be determined by 
the change of a single vortex solution in a widely sepa- 
rated vortex-lattice configuration with a lattice spacing 
dh [E 0- The variation of the single vortex solution 
occurs dominantly around the boundary of two vortices 
~ di/2, proven to coincide with an asymptotic solution 
of the single vortex. In this respect we proceed as follows. 
First, we find the asymptotic solution of a single vortex 
away from the vortex core. Second, we show that the 
variation of the vortex solution is given by the asymp- 
totic solution. Third, we represent the vortex interaction 
in terms of this solution. 

We introduce the following ansatz for an asymptotic 
solution of the single vortex configuration 

s(x, z) = rj v (r, z) - i] s (z) = A(z) R(r) , 

q t (x, z) = Q\ (r, z) - Q\(z) = B(z) R(r) , 

q t (x, z) = Q\(r, z) - Q\(z) = C{z) K^r) , (4) 

where the superscript v represents the single vortex solu- 
tion and s denotes the uniform solution with the radial 
coordinate r or rectangular coordinates Xi in two dimen- 
sion. Then, Eq. <|3j> becomes 

(V 2 - k\) R(x) = 0, (V 2 - k\) K t {x) = 0, 

/' 2x 2 Qf\ Al 27fQ? 
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(5) 



in q z = gauge. It is straightforward to see R(r) = 
Koi^ir), lZe(r) = Ki(n2r), and 1Z r {r) — in polar coor- 
dinates. Here, K\ and K2 are constants for separation of 
equations. Scaling the radial coordinate by f = K%r, we 
find that Eq. ([5]) can be rewritten in terms of only a single 
parameter n = For other equations, we need to 

solve them numerically, taking the regularity conditions 
at the horizon. It turns out that resulting solutions de- 
pend on k and a, b, c, which are defined at the horizon, 
a = A(l), b = £>'(!), and c = C(I). In addition, we find 
that such solutions are characterized only by bj a and k 
due to the scaling symmetry of Eq. ([5]) . See appendix A 
for the numerical analysis [22j. 

Having the asymptotic solution, we evaluate the ef- 
fective interaction between vortices in the dilute vortex- 



lattice configuration p], 0| • We introduce 

r, = rf + 5t] , Q t = Q v t + SQ t . 
Q l = Q v l + 5Q A , = n6 + 6<f> 



(6) 



where the solution with the superscript v represents the 
single vortex configuration in a Wigner-Seitz cell, while 
the "<5" part expresses the variation of the single vor- 
tex configuration around the boundary of the Wigner- 
Seitz cell. n is the winding number of the vortex. 
6<f> = J2i^o n arg(a? — afj) is chosen for a multi- vortex con- 
figuration, where x^ is the core position of each vortex. 
Srj and SQ t would be much smaller than rf and Q v t inside 
the Wigner-Seitz cell, respectively. On the other hand, 
it is not obvious if SQi is much smaller than Q\ near 
the boundary of the Wigner-Seitz cell because Q1 will be 
also small. However, it is natural to expect that Sri and 
SQt are much larger than SQ 2 near the boundary As 
a result, we obtain the following linearized equations of 
motion near the boundary 



+2^SQ t - 2^5Q t = 0, 
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Sr] = 0, 

"+ 7 5 Q'i + (- - ^Q* - = o. (7) 

An important aspect is that these equations are essen- 
tially the same with those for the asymptotic configu- 
ration of the single vortex, valid when rf 3> Srj and 
Qt > SQt with rf or Q v t > SQ 2 . This property leads us 
to write down the variation of the solution in terms of the 
asymptotic solution for the single vortex configuration 



Sr] ■ 
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{x-x h z). (8) 



Expanding the action ([TJ around a vortex solution to 
second order and using equations of motion ([3]), we arrive 
at 



Sil {2) = i [dzcPxd^ 



difP l! (r) v + -Srj 



5Q„( 
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-SB^ 1 
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We observe that only surface terms contribute to the cor- 
rection for the grand potential, where these boundaries 
correspond to the AdS4 boundary [z = 0), the horizon 
(z = I), and the boundary of the Wigner-Seitz cell. The 
regularity condition on the horizon does not allow con- 
tributions from the horizon. In addition, the dilute vor- 
tex configuration guarantees that the contribution along 
the AdS4 boundary is much smaller than that from the 
boundary of the Wigner-Seitz cell [22j |. Therefore, the 



FIG. 1: Asymptotic vortex solutions for dimension 2 conden- 
sation, charge density, and magnetic flux, respectively, where 
k = 1 is used. Both the blue line with [a/b = — 1, a/c = 2] 
and the green line with [a/b — 1, a/c — 6] correspond to re- 
pulsive interactions between vortices while the dashed red line 
with [a/b — — 1, a/c — 6] results in attractive interactions. 



FIG. 2: Effective interactions between vortices as a function of 
k. With decreasing n the effective potential changes from re- 
pulsive to attractive. The nearest interaction was considered 
where d L = 3. Left: T/T c = 0.907, b/a = -1, a/c = 3 (Blue 
Solid), and a/c = 10 (Red Dashed). Right: T/T c = 0.388, 
b/a = —1, a/c = 3 (Blue Solid), and a/c = 10 (Red Dashed). 



relevant contribution is from the correction at the bound- 
ary of the Wigner-Seitz cell, which is attainable by the 
asymptotic solution (|8|) . Finally, we obtain the change of 
the grand potential for the i th cell, 

5tlW~- f dz Idl h- \irSriV(ri v + -6r}) 
T Jo Ji ^z 2 2 

' SQ t V(Q v t + ^SQt) + SQxVx(Q v + ^SQ)} 



m 

= irE [ CR o(n) -(.A- B)K Q {V2Kr-S\ (10) 

with A =J^dz^-, B =f£dz^£- and C =J^dzC(z) 2 . 
For the analytic expression in the last line of Eq. (flOl) . 
we used an identity in Ref. [23j . 

It is possible to understand the physical meaning of 
Eq. (fTU|) . The interaction potential consists of both first 
order and second order contributions in "<5" , where the 
former represents interactions between the i th = vortex 
and others i th ^ 0, and the latter expresses those between 
other vortices i th ^ except for the i th = vortex. 

This expression is formally identical to the effective 
interaction between vortices in the Abelian-Higgs model, 
where the first term results from the variation of the su- 
percurrent while the second originates from that of the 
Higgs field around the boundary [l], Q ■ An important in- 
gredient is that coefficients of the vortex interaction are 
given by integrals in the z-direction. In addition, the 
k dependence of the interaction potential is much more 
complicated since such coefficients are functions of the 
parameter k. In this respect the role of the parameter k 
is not completely clear yet although tuning k results in 
the change of the vortex interaction. 

Figure [1] shows dimension 2 condensation, charge den- 
sity, and magnetic flux for the asymptotic single-vortex 
configuration, respectively. It is interesting to observe 
that when U(l) charge density decreases rapidly near the 
vortex core, the effective interaction between vortices be- 
comes more repulsive. As long as b/a remains positive, 
we do not see any change from repulsive to attractive 
interactions. In this case the system might lie in a deep 
type II regime. Therefore, we focus on b/a < hereafter 
to study the superconducting transition. 

We classify our systems into four classes under the 



condition of b/a < 0, depending on the density of U(l) 
charges and the ratio of a/c. First, we fix the density 
of U(l) charges, determining the chemical potential. We 
expect that the regime with a/c <C 1 belongs to the type 
II superconductivity because the first term of the vortex 
interaction in Eq. (|10l) becomes larger than the second 
term, resulting in repulsive interactions. Physically, this 
relation implies strong supercurrents around the vortex 
core, consistent with the picture of type II. On the other 
hand, the regime with a/c » 1 will show that interac- 
tions between vortices change from repulsive when k> 1 
to attractive when n <C 1. Notice that k is introduced 
into the second term, reducing it with k> 1 and enhanc- 
ing it with k<1. Both A and B are positive definite, 
decreasing monotonically as we increase k. We uncover 
that the regime with b/a < gives rise to A — B > 0, 
allowing the possibility for the change of interactions. 
Fig- El confirms our expectation, that is, interactions be- 
tween vortices become attractive when n < n t , where ttt 
can be regarded as the tricritical point. 

Next, we consider cases with a fixed k. When k is 
rather large, it is difficult to find the tricritical point /i t , 
originating from smallness of the second term. In this 
respect it is better to start from a small enough k. Then, 
the effective interaction is attractive when /i > fi t while 
it becomes repulsive when [i < nt ■ Figure [3] shows a sur- 
face of tricritical points in the space of (k,/x/T) with a 
fixed T, a/c > 1, and b/a < 0, where effective interac- 
tions between vortices vanish exactly. The vortex inter- 
action is attractive inside the ellipse while it is repulsive 
outside the ellipse. We claim that this ellipse serves a 
general criterion for the fluctuation-driven first-order su- 
perconducting transition in strongly coupled conformal 
field theories, possibly occurring in the vicinity of quan- 
tum criticality. 

In this study we try to answer how to classify strongly 
interacting field theories, considering the nature of the 
superconducting transition. The holographic supercon- 
ductor model is our main ansatz as an effective low energy 
theory, expected to describe certain classes of strongly 
coupled conformal field theories. The effective interac- 
tion between vortices is our central object, allowing us 
to distinguish the type II superconductor from type I, 
where the former will show the second order transition 
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FIG. 3: Tricritical surfaces in (b/a,a/c,K) with a fixed [a, for 
the left panel and (/j,,a/c, k) with a fixed b/a for the right 
panel, respectively, where vortices do not interact with each 
other. The interaction potential is attractive inside the ellipse 
while repulsive outside it. 

while the latter will display the first order. As shown, 
an asymptotic solution for a single vortex configuration 
plays an essential role for the effective interaction. The 
effective interaction between vortices turns out to be a 
complicated function of both k and /i/T, where the pa- 
rameter k is introduced to play basically the same role 
as the Ginzburg-Landau parameter. We find a surface 
of tricritical points in the parameter space of (k, /z/T), 
where the effective interaction vanishes, which separates 
the first order from the second order, proposed to be a 



general criterion in classifying quantum critical metals. 

There are various unsolved questions in this direction. 
First of all, a possible topological term such as the ax- 
ion term |24j may play an important role in the vortex 
interaction. It can assign the U(l) charge to a vortex, 
modifying their interactions. We suspect the possibility 
of the BKT transition [25[ , resulting from their Coulomb 
interactions due to the assigned U(l) charge, where 1/q 2 
in the momentum space becomes In r in two space dimen- 
sions. In addition to this problem, the role of the pairing 
symmetry is not investigated, where non s-wave super- 
conductivity arises in strongly interacting electrons (tJ- 
Furthermore, it should be studied the role of fermions in 
the vortex interaction. 
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Appendix A: Numerical analysis for Eq. (5) 

We start with discussion about a uniform solution in the original holographic superconductor model [S. A. Hartnoll, 
C. P. Herzog, G. T. Horowitz, Phys. Rev. Lett. 101, 031601(2008)], where Q s t (z) and r) s {z) depend only on the 
z-coordinate. It is straightforward to derive equations of motion from Eq. (3) 



•/' + (H>' + (w + i F> = - *'-57*-* (A1) 

The regularity at horizon (z = 1) gives the following conditions 

Q t (l) = 0, Q' t (l)=v, Q»(1) = _Y, (A2) 

., . 2 ... . 1 / 8 v 2 
7]{1)=U, 7]{1) = -U, J] (1) = -u I -- - — 



(A3) 



Using the above conditions, one can find solutions in terms of u and v, based on the shooting method. Near the 
boundary [z = 0), the solutions behave such as 

i](z) ~ Oi(u,v)z + 2 {u,v)z 2 ^ , Qt{z) — fi(u, v) + p(u, v)z H . (A4) 

When we are considering operators with dimension 1 or 2, we should constrain solutions with either O2 = or 0\ = 0, 
respectively. Therefore u and v are not independent but related with each other. As a result, the space of solutions 
becomes one dimensional, allowing us to take 'V as a parameter for the solution of the holographic superconducting 
state. In other words, v controls cither temperature or charge density of the system. According to the AdS/CFT 
dictionary, the total charge is given by Q <~ — J d 2 xaQ' t (0). When T = |^ is fixed, the charge density varies as a 
function of v. If the charge density is fixed, temperature changes as a function of v. One can say a similar statement 
for the chemical potential, a Qt{0). Inserting the uniform solution of Qf (z) and rj s (z) into Eq. (5), we can solve them 
numerically. 

Equation (5) has two parameters of K\ and Ki- Performing the scaling as discussed in the manuscript, we obtain 
the following regularity conditions near the horizon for A(z), B(z) and C(z), 

/2 2k 2 \ 4a buv av 2 Aan 2 2aK 4 . , „. 

A(l) = a, A'(l)=a(- + — J, A»(l) = --- — - — + — + —, (AS) 

B(1) = 0, B'(l) = b, B"(l) = -\auv + \b(-u 2 + 2n 2 ) , (A6) 

0(1) -C C(l)-c(i-^). C- ( l) + £ + (AT) 

It is straightforward to see the scaling symmetries in Eq. (5). Equation for C(z) remains invariant after scaling as 
cC(z) with a parameter c. In this case c = 1 is allowed due to the boundary condition in Eq. (A7). Equations for A(z) 
and B{z) also allow scaling, unchanged after aA(z) and aB(z). Therefore b/a and k arc only relevant parameters, 
governing equations for A(z) and B(z). As a result one may regard the asymptotic solution of the single vortex 
configuration as 

s = aA(v, k, -, z\Kq(_\/2kT), q t = aB(v,K,,-,z)K (V2~Kf), qe = cC(v, z)K 1 (f) , (A8) 
v a / v a / 

where f means a re-scaled coordinate with K2. We emphasize arguments in each function. 
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Appendix B: Derivation of the variation for the grand potential 



Inserting Eq. (6) into Eq. (3), we find the following linearized equations 

V 2 8r) - m 2 6f] ~ Sj]Q v2 - 2rj"Ql8Q> 1 = 0, 



(Bl) 



V^B^ 



2rf5i 1 Q v - r\ 



v2. 



(B2) 



proven to be valid near the center of a vortex. One can see that this approximation is reasonable only when Srj and 
SQt are both larger than SQ 2 and smaller than rf and Q\ . The boundary of a Wigner-Seitz cell also satisfies these 
conditions. In this respect the linearized equations are valid not only near a vortex but also the boundary of the cell. 
This is a simple extension of the observation in L. Kramer, Phys. Rev. B 3, 3821 (1971). 

Inserting the vortex solution [Eq. (6) with Eq. (8)] into the effective gravity action [Eq. (1)] and expanding the 
action to the second order, we obtain the following expression for the change of the grand potential in a cell 



SCI = i / dzd 2 xd IA 



tiifD" | if + -Sr) 



B v » v + -5B^ V 
2 



(B3) 



where Eqs. (Bl) and (B2) are utilized. In this derivation we need to worry about singular parts from X^ v . 8X^ V 
vanishes identically because the singularity appears completely outside the cell i. The only term that we have to 
concern is 8B IIV X ,1V , however this turns out to vanish when we are considering the configuration of SB^ = at the 
origin of the cell. 

Changing Eq. (|B3[) into surface integrals, we have three kinds of boundaries. The first is the boundary at the 
horizon (z = 1) of the black hole and the second is that of the AdS space (z = 0). The last is the boundary of the 
Wigner-Seitz cell. The first contribution vanishes identically thanks to regularity conditions at the horizon. For the 
z = boundary, the contribution must be considered carefully. Actually, this contribution could be important, when 
a distance between vortices is comparable to a size of a vortex. However, we are taking the dilute gas limit, thus the 
variation from the single vortex solution will be concentrated on boundaries of Wigner-Seitz cells. 

The surface integral for z = is given as follows 



d 2 x y 



Sr]V z (y 



B vzv + \^B ZV 



(B4) 



In the dilute limit Srj and SQ^ have nonzero values only near the boundary of a cell. Thus, the integration range is 
effectively small. As positions of vortices are far from each other, this contribution almost vanishes and it is much 
smaller than the third contribution given by the integration along the z direction at the boundary of Wigner-Seitz 
cells. This dilute approximation serves the validity of our calculation. Therefore, our correction of the grand potential 
is well approximated as 



sn 



i J dz f dlfiiy/=g 5 V V 1 (rf + ^r?) + 5Q V (b v iv + X -8B l 



(B5) 



where hi is a unit vector orthogonal to the boundary of the Wigner-Seitz cell. This leads to Eq. (9). 



Appendix C: Derivation of Eq. (9) 



In this section we will derive the following formulae 

Y,<fdlh-K (\r-r l \)e(r-r i )xWx (^K (r)9 + ± £ 7f (|r - r j\)6(r- f,)) = Tr^ifoG^I), (CI) 

Y,fdlh-K {\r-r l \)v{K {r)+ l -Y, K ^\ r ~ r i I)) = E WiD- (C2) 



i#0 ' 
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For convenience, we define a fictitious coordinate w and vector fields, hi(r) = K$(\r — ri\)w, where is a center of a 
lattice. Then, the above equations can be written as follows 



Fl = J2f d S 



i#0 



x ( V x hA > = ir^hi 



■ w, 



F2 



dS ■ hi x V x [h 



(C3) 
(C4) 



where dS* is an area element orthogonal to the boundary surface of a cell, i.e, dS* = dl dw n. 
Using the divergence theorem and taking integration by parts, one can rearrange Fl into 

fi = 1 d3x { v x { h ° + \ Yl ^)'( v x h ^ - { h ° + /l j)'( V x V x fti) i 



Y / d M v x r° + 2 



(C5) 



where we have used V X V X hi + hi = 2ir5 2 (r — fi)w with some algebra. Actually, the second term in the last is equal 
to —Fl. In order to show this, we consider the following combination, 



2 ^ ^ 



Y I dS ■ < (h + ~ ^ hj^j xV x hA + Y i dS ■ IhiXV x (h + ~ ft, 
= i" • V J2 h o ■ ( h o + ~ E h 3') 

= IdS-wlh k -h + (h + h k )- hi + \j2 hi ' Y ft i ' 



(C6) 



where fc means one of other vortices. The integrand is symmetric under interchange of ho and h k . Now, we may 
take rfe as one of the nearest neighbor. Then, the symmetry means that the integrand is a vector whose direction 
is along the boundary surface of the cell. Thus, the above combination should vanish, leading us to conclude that 
Fl = tt^^o ' This argument is in parallel with that in L. Kramer, Phys. Rev. B 3, 3821 (1971). For F2, 
one can see that F2 = —Fl from Eq. (|C6[) . This completes our derivation. 

I 



